Spatial dispersion in CaF2 caused by the vicinity of an excitonic bound state. 
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The microscopic mechanism beyond the optical anisotropy of an ionic crystal which occurs for 
short wavelengths is investigated. The electron-hole, two particle propagator and its analytical 
behaviour close to the band edge of the one particle continuum plays a major role for the mechanism 
of this optical anisotropy. Especially for an ionic crystal the two particle bound state, the exciton, 
is of special importance. In this way we argue that the so called "intrinsic birefringence" in CaF2 
is neither intrinsic to the material nor it is birefringence. Instead it is spatial dispersion caused by 
the vicinity of a dispersive optical absorption given by the excitonic bound state. We propose a 
model which connects the bound state dispersion with the band structure and a model potential for 
a screened coulomb interaction. Based on these considerations we predict a wavelength dependence 
of the dielectric function approaching close to the bound state level e ~ (A — Ao)~^, where Ao is the 
wavelength of the excitonic bound state level. 

PACS numbers: 



I. INTRODUCTION 

There is no doubt that CaF2 is the key component 
for estabUshing 157nni UV lithography. Its perfect cu- 
bic symmetry together with a great homogeneity of large 
single crystals which are very stable under normal condi- 
tions makes it a perfect material for homogeneous optical 
lens design. 

However, there is a severe shortcoming of the mate- 
rial on which the attention was focused within the. last 
year. This is the so called "intrinsic birefringence" .El For 
the first glance it is astonishing. Looking at standard 
textbooks (e.^) we note that a cubic crystal does not 
show any birefringence. On the other hand experiments 
observe a spatial dependence of the dielectric properties. 
Light propagating along different directions through the 
crystal shows a deviation between different polarization 
directions which is of the order of 10~^ nm/cm for UV 
light with 157nm. Or with other words the dielectric 
function depends on the direction in which the light prop- 
agates through the crystal and shows a difference which 
is of the order of 10~^. 

How can this obvious misfit between standard text- 
books and recent observations in CaF2 be resolved? The 
answer leads back to a ph;j£sical phenomenon which is 
known as spatial dispersionjj 

The statement "a cubic system does never show bire- 
fringence" can be proven exactly, since in a cubic system 
the ellipsoid corresponding to the dielectric function is 
always a sphere. The question remains: What is a "sys- 
tem"? The crystal of CaF2 has exactly cubic symmetry 
but the "system" of light interacting with this perfectly 
cubic crystal shows a deviation from cubic symmetry. 
This is due to the fact that the light carries a q vec- 
tor pointing in the vacuum in the direction of propa- 
gation. This is a consideration Wjiiich has been already 
made more than one century agoQ. However the speed 
of light is much larger than any speed in a solid state 
system and therefore one can usually safely neglect the 



wave- vector dependence of the propagating light. This 
statement is identical to the wavelength of light being in 
the optical window much larger than any length scale in 
a solid. 

There are however two important exceptions where the 
assumption made above ceases to be valid. The first (I) 
is when one is by far leaving the optical window. This 
is the case at X-ray scattering experiments. Here the 
wavelength of light is of the order of the lattice constant 
and this causes a wave- vector (spatial) dependence of the 
resulting scattering signal. The second case (H) which is 
valid for CaF2 and many other solids is the case of light 
in the vicinity of a dispersive absorption occurring also 
at small wavevectors (large length scales). Since the di- 
electric function is an analytic function it has to follow 
causality which is expressed in the Kramers-Kronig rela- 
tion between its real and imaginary part. The vicinity 
of a band edge can cause spatial dispersion in a semicon- 
ductor system (see e.g.El). 

In a semiconductor, the absorption is mainly domi- 
nated by the band edge and an exciton, the two particle 
bound state of the electron-hole pair only plays a minor 
role. The situation changes if we look at ionic crystals 
like CaF2 . Here screening of the bare coulomb interaction 
is less effective which expresses itself in deep excitonic 
bound states. This is expressed in the exciton intensity 
being inverse proportional ta,the third power of the static 



dielectric function / 



The connection between a 



dispersive excitonic bound state and spatial dispersion 
has been already discussed for KLQ 

When discussing possible alternative material designs 
to CaF2 it is important to understand the microscopic 
mechanism beyond the spatial dispersion in great detail 
and it is in particular important to know and predict the 
exciton dispersion. 

The paper is organized as follows: In section II we 
define the system in the case when the approximation 
of locality breaks down. In section [II we review the 



connection between the electronic part of the dielectric 



a) 



function and the bare two-particle propagator. We look 
at the free two particle propagator in the absence of any 
screened coulomb attraction and its connection to the 
band structure. The connection between the band struc- 
ture, which can e.g. be based on DFT (density functional 
theory) calculations, and the dynamic polarization is in- 
vestigated in section IV . In the following (section ^ we 
discuss analyticity based on the resulting pole structure 
of the dielectric function. This describes to a large part 
the situation in a semiconductor where no deep bound 
state is present. Section ( |Vl[ ) serves to introduce the two 
body problem of the electron hole pair which leads to the 
excitonic bound state when a reasonable assumption for 
the weakly screened coulomb attraction has been made. 
Based on the model assumptions made, we discuss the 
spatial dis pers ion of CaF2 close to the bound state level 
in section VII and compare our predictions with the ex- 
perimental work of.EJ 



II. SYSTEM OF CAF2 INTERACTING WITH 
LIGHT 
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Light has a q vector pointing along its propagation di- 
rection. This wave-vector becomes important when the 
light is interacting with a long wavelength (q — > 0) exci- 
tation of the system. In an isotropic medium this wave 
vector dependence will only depend on the absolute value 
of q. In liquids and glasses this leads e.g. to the obser- 
vation of hydrodynamic. modes (phonon modes) in light 
scattering experiments. eI In single crystals the situation 
is different as not only the dependence on the absolute 
value of the wavevector of the light is important but also 
the direction of the wave vector with respect to the crys- 
tal orientations. The crystal of CaF2 possesses ideal cu- 
bic symmetry. That means in the static limit that also 
all response functions have to show this symmetry. Dy- 
namically there are situations especially close to a strong 
excitation at long wavelengths, when the breaking of cu- 
bic symmetry due to the incident light gets important. 
Therefore not only the dynamic of a dielectric response 
has to be considered but also the wave-vector depen- 
dence. The electric field E is connected with the electric 
induction D in a non-local way. 

D(q,c^) = e(q,c^)E(q,c^) (1) 

where e(q, oj) is the dynamic dielectric function. For the 
physics of the spatial dispersion in CaF2 we will focus 
on the part of the dynamic dielectric function which is 
caused by electronic excitations. In general the dielectric 
function is an analytic function defined on the complex 
cj-plane. Along the real axis it has an real and imaginary 
part. 

e(q,w) = e'(q,c^) + ze"(q,c^) (2) 

Only in the static limit uj — > one can assign a physical 
value to the real and imaginary part separately. The real 



FIG. 1: Schematic plot of the pole structure of the two 
particle correlation function n*'(q,tj) as it can occur for dif- 
ferent directions of the wave vector. n''(q, oj) is in the case 
when no excitonic bound state exists proportional to the di- 
electric function. In Fig. a) we have plotted the band gap. 
The band gap can have a different shape when looking along 
different q directions. Above the band gap the two particle 
continuum starts where there is always a non zero imaginary 
part of the dielectric function. The light dispersion uj = cq 
is plotted schematically. In Fig. b) the real and imaginary 
part of the dielectric function resulting from Kramers-Kronig 
relations are plotted in a schematic way. When approaching 
the band gap the real part of the dielectric function bends up 
(but will not diverge in a 3D system) and the imaginary part 
will show a yjui — uio behaviour. 



part is connected with polarizabilities and the imaginary 
part with mobile charges (e"(aj) = limi^^o 47rcr(w)/a;) 
which are absent in an insulator. For a general crystal 
the dielectric function is a tensor. Since we want to keep 
the equations as transparent as possible we do not write 
tensor indices. We assume a coordinate system placed 
along the main axes of the cubic crystal. For a cubic 
crystal like CaF2 we can concentrate on specific direc- 
tions of the q vector. Symmetry considerations allow to 
decompose all tensorial quantities into its symmetry com- 
ponents. This will provide rules to fill all the components 
of the tensorial dielectric function. 
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III. THE DYNAMIC DIELECTRIC FUNCTION 

The optical frequency regime is above the infrared op- 
tical phonon modes and below excitations over the band 
gap. In this regime the dynamical dielectric function can 
be decomposed in two parts. 



e(q,w) 



h+ 



(3) 



Where Coo is the IR-dielectric function which in a cubic 
crystal like CaF2 will be not only real but completely 
isotropic. The subscript c» refers to frequencies above 
the IR optical phonon modes but is still far below exci- 
atations over the bandgap. Wave- vector dependence only 
enters from the second part which stems from electron 
hole excitations. The electron-hole dielectric function 
is proportional to the proper electron-hole polarization 
propagatorO 



(q,c^)^n*(q,c^) 



(4) 



Where 11* (q, uj) is the propagator of an electron hole pair 
with total momentum q. In the following we do not want 
to focus on details but rather want to outline the pole 
structure of this two particle propagator since it gives 
further insight into the problem. 



IV. CONNECTION TO THE BAND 
STRUCTURE 

If we look for the first glance at the bare electron hole 
propagator n*'(q, w) we can make the connection to the 
band structure. 



k ^ 



/(k, q,w) 



-(k + q)-i?''+(k)) 



(5) 



The detailed form of the numerator can be taken from 
textbooksa while the pole structure is connected with the 
electronic band structure. In general Eq. (^) should also 
contain a sum over all bands. Since we are approaching 
the band gap from the low frequency site and we know 
that CaF2 has a direct band gap we are for a first glance 
only interested in the highest valence band (k) and 
the lowest conduction band (k) around the F point 
at k = (0,0,0). Even in a cubic crystal the band struc- 
ture will show a diff'erent k dependence when looking 
along different directions in the crystal. In particular the 
(1, 0, 0) direction wiU differ from the (1, 1, 0) and (1, 1, 1) 
direction. 



V. ANALYTIC PROPERTIES OF THE 
DIELECTRIC FUNCTION 

From the analyticity of II" (q, oj) we can derive the 
qualitative behaviour of its real and imaginary part. This 




FIG. 2: Schematic plot of the pole structure propper po- 
larization n*(q,t<;) (proportional to e^- (q, w)) as it can 
occur for different directions of the wave vector in an ionic 
crystal where a strong excitonic bound state is present. In 
Fig. a) we have plotted the band gap. The band gap can 
have a different shape when looking along different q direc- 
tions. Above the band gap the two particle continuum starts 
where there is always a non zero imaginary part of the dielec- 
tric function. The bound state can occur below the one parti- 
cle continuum. For simplicity we have only plotted the ground 
state level of only one exciton. The light dispersion ui = cq is 
plotted schematically. In Fig. b) the real and imaginary part 
of the dielectric function resulting from Kramers-Kronig rela- 
tions are plotted in a schematic way. When approaching the 
bound state the real part of the dielectric function diverges 
as {lo — ujo)~^ where ujq is the bound state level. 



is plotted schematically in Fig. |l|. In Fig. |i|a the struc- 
ture is shown in the uj — q plane. There is the large band 
gap. An excitation over the bandgap is the creation of 
an electron-hole pair. Above the band gap there is a 
continuum of one particle excitations. The light disper- 
sion is shown as well. Since the speed of light is much 
larger than any electronic speed one can usually in the 
optical frequency range safely neglect the wave- vector de- 
pendence of the light and use q = 0. This breaks down 
when excitations are considered which come close to the 
band gap. In this case it can become important that the 
band curvature is different along different directions in 
the crystal. Schematically we have plotted here a ficti- 
cious (1, 1,0) and (1,0,0) band curvature. 

In Fig. |l|b the resulting real and imaginary part of 
n'^(q, cj) is plotted resulting from the Kramers-Kronig 
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relations. In a three dimensional system the imaginary 
part of the one particle continuum shows a y/uF^^ojQ be- 
haviour. Therefore the real part does not diverge at the 
band edge (opposite to a ID or 2D system) but increases 
up to a finite value. In order to keep the plots in Fig. |l| 
as simple as possible we have not plotted the real part of 
the polarizability above the gap. 

If for different wave- vector directions in the crystal the 
band gap is shifted one can observe a spatial dispersion. 
This is the case in semiconductors where no deep exci- 
tonic bound state existsJj Here the spatial dispersion is 
caused by the wave- vector dependence of the one particle 
continuum. A weak attractive electron-hole interaction 
will not change the picture as long as no strong excitonic 
bound states are created below the one particle contin- 
uum. Therefore in these cases the dielectric function can 
already be described using e ^ Hq. The situation is com- 
pletely different in a strongly ionic crystal like CaF2 here 
the excitonic bound state becomes important. 



n°. This crossing marks the bound state level where 
a new pole occurs. Such a situation with a bound state 
is schematically sketched in Fig. ^j. Below the one par- 
ticle continuum the exciton level occurs as a true pole. 
The dynamic polarization and therefore also the dielec- 
tric function has to follow Kramers-Kronig relations con- 
necting its real and imaginary part. In the vicinity of a 
true pole the real part of the dielectric function has to 
follow the asymptotic behaviour as 



~h+ 



(q,c.)~5R(n*(q,c.)) 



1 



27rc 



=(q) 



.(q) 
(8) 



The wave- vector dependence of the exciton level a;ea:(q) 
(or Aea;(q)) wiU in general differ along different directions 
in the crystal. This is schematically plotted in Fig. ||a. 
Exactly this difference causes the spatial dispersion seen 
in ionic crystals like CaF2 and BaF2. 



VI. THE TWO PARTICLE BOUND STATE OF 
THE EXCITON 

In a strongly ionic crystal screening is less effective and 
the coulomb attraction between the positively charged 
hole and the negatively charged electron leads to the for- 
mation of a strong two particle bound state. The solution 
of this two body problem involves in geaeral a solution 
of the Bethe-Salpeter equation (see e.g.E2l). However the 
basic principal of this bound state formation can already 
be understood with a much simpler interaction. The full 
solution for a complicated electron Jiple interaction for 
CaF2 has been taken into account intll here we just want 
to highlight the basic principles under which such a solu- 
tion leads to an excitonic bound state. When we assume 
the attractive interaction to be only a function of the to- 
tal momentum of the electron-hole pair C/(q), the equa- 
tions can be solved easily and one gets a RPA (random 
phase approximation) denominator for the proper polar- 
ization which stems from a summation of a geometrical 
series. 



H*(q,c^) 



n°(q,c.) 



l-C/(q) HO(q,c^) 



(6) 



This is also known as the ladder approximation to the 
Bethe-Salpeter equation. New physics will occur when 
the denominator of Eq. (||) shows additional poles beside 
the one particle continuum. This is the case when the 
denominator of Eq. (||) gets zero. 



l-C/(q)n"(q,c^) = 



(7) 



If this happens outside the one particle continuum a true 
bound state, the exciton is formed. Finding the bound 
state level can be schematically seen in Fig. |l|b. When 
drawing a hne at /(w) = 1/U{q) it can for large t/, 
for large attractive interactions cross the real part of 



VII. SPATIAL DISPERSION IN CAFa AND 
PREDICTIONS FOR EXPERIMENTAL 
OBSERVATIONS 

In pview of the measurements performed by Burnett 
et alJ^I we need the refractive index along different di- 
rections. In an isotropic medium this is where the di- 
electric function is diagonal with three identical diagonal 
elements this would be simply given by: 



n{Lo) = e{uj) 



(9) 



with e being the (real part of the) diagonal part of the 
diagonalized dielectric function. 

The concept of a refractive index becomes by far more 
sophisticated as soon as the medium shows anisotropies. 
This is also the case for CaF2. Here one has to go into 
a new coordinate system whose z-direction is determined 
by the q-direction of the incident light. Two possible 
directions perpendicular to the z-direction complete this 
new coordinate system and are denoted with Greek let- 
ters a, (3. Within the a-/3-plane lies the polarization vec- 
tor e pointing in the direction of the electric induction 
D of the light. Within this new coordinate system the 
calculation of the refractive index is given by the zero of 
a two dimensional determinant .u 



Det 



'^5o,,p - e^j3(q,w) 







(10) 



In general the eigenvalue problem of Eq. ( pi]| ) has two 
solutions which correspond to the main polarization axes. 

There is a further difference to an anisotropic medium. 
This is the fact that the wave-vector q is not in general 
parallel to the Poynting-vector S which points along the 
direction of the group velocity which is the analogon of 
the path in classical geometrical optic. This is expected 
to be small in a cubic crystal. 



4 



CaF2 intrinsic birefringence 
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In Fig. ^ we have plotted the wavelength dependence of 
An as it would result from the esfjeriiientally observed 
bound state level lying at 112nm£^ 

It can be seen that describing the experimental obser- 
vations with Eq. ( p^ (dotted line) giv3 
ment with the experimental observati: 
just a 1/A^ bfeMTOtJ-^traasIhgd ii|rc) as it would be ex- 

ftf^)^'^(l-F^'i^^i^^'ileT^^'!§^'?^^"— dielectric func- 
tion. WcTthcr'cforc clairri in this work th at the anisotropy 



(q) 



(12) 
(13) 



s a better agree- 
m (crosses) than 



FIG. 3: Comparison of the data fromcl (crosses) with a 1/A^ 
(dashed line) and a 1/(A — Aea;(q)) (dotted line) behaviour. 
The behaviour results from a Taylor expansion of the 

dielectric function (dashed line). A comparison with the func- 
tion which results from the vicinity of a dispersive bound state 
(exciton) is plotted (dotted line). iFsr the exciton wavelength 
the 112nm known from literatural3 were used. The exciton 
model produces a better description of the data at smaller 
wavelength. It further makes the prediction that the dielec- 
tric properties will for wavelength below 157nm diverge as 
1/(A — Aea;(q)) when aproaching the excitonic bound state. 



So the refractive index is dependent on q and can for 
a given experimental setup be dependent on the polar- 
ization direction of the incident light. 



n — n(q, oj) — n(q, 27rc/A) 



(11) 



In a cubic crystal like CaF2 there are several exceptions 
for the wave-vector q and the Poynting vector pointing 
along the same direction where no dependence on the po- 
larization direction occurs as one can derive from symme- 
try considerations. One is for example light propagating 
along the main axis of the crystal (e.g. q = (1, 0, 0) now 
measured in the coordinate system of the crystal axes) . 

The dielectric function diverges close to the excitonic 
bound state depending on the wave-vector q according 
to Eq. (||) with Aea;(q) is the wavelength of the exci- 
ton for the given q vector. This wavelength must be 
for CaF2 closeJx) the 112nm where the exciton is found 
exp er iment allyllj . 

Note that Eq. (||) is quite different to what would re- 
sult from just a formal Taylor expansion of the dielectric 



which S^Seen in Cal^Oat low wavS69igths is caused 
nrrljy the exciton dispersion of an excitonic level around 
112nm (11.2eV). Especially measurements in the wave- 
length range between 112nm and 157nm we predict to 
show a clear 1/(A — Ae:r(q)) behaviour. 



VIII. CONCLUSION 

We argue that the optical anisotropy observed in CaF2 
is caused by exciton dispersion. The delocalized exciton 
lies at approximately 112nm. Using this excitonic model 
the measured data on the optical anisotropy can be well 
explained and it is predicted that the dielectric function 
should diverge for lower wavelengths as 1/(A — Aea;(q)). 
Here Aea;(q) is the wavelength of the excitonic bound 
state. 

If the connection to the exciton is indeed the domi- 
nant mechanism beyond the optical anisotropy all possi- 
ble compensation approaches have to be reviewed. E.g. 
forming solid solution crystals, if possible at all, would 
only lead to an improvement, if the exciton dispersion 
will be smaller in such a crystal. Further the subtle dif- 
ference between the spatial dispersion and birefringence 
has to be taken into account in a very careful way when 
performing lens design. 
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